The neutron electric dipole form factor in the 
perturbative chiral quark model 



Claudio Dib|, Amand Faessler$, Thomas Gutsche$, 
Sergey Kovalenkof , Jan KuckeiJ, Valery E. LyubovitskijJ, 
Kem Pumsa-ard| 

f Departamento de Fisica, Universidad Tecnica Federico Santa Maria, 
Casilla 110-V, Valparaiso, Chile 

| Institut fur Theoretische Physik, Universitat Tiibingen, 
Auf der Morgenstelle 14, D-72076 Tubingen, Germany 

E-mail: cdibOf is .utf sm. cl , amand.faessler@uni-tuebingen.de, 
thomas . gut scheOuni-tuebingen . de , Sergey . Kovalenko@usm . cl , 
kuckeiOtphys . physik . uni-tuebingen . de , 
valeri . lyubovit ski j Ouni-tuebingen . de , 
pumsaOtphys . physik . uni-tuebingen . de 

Abstract. We calculate the electric dipole form factor of the neutron in a pertur- 
bative chiral quark model, parameterizing CP-violation of generic origin by means 
of effective electric dipole moments of the constituent quarks and their CP-violating 
couplings to the chiral fields. We discuss the relation of these effective parameters to 
more fundamental ones such as the intrinsic electric and chromoelectric dipole moments 
of quarks and the Weinberg parameter. From the existing experimental upper limits 
on the neutron EDM we derive constraints on these CP-violating parameters. 
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1. Introduction 

The origin and manifestations of CP violation have been important issues in particle 
physics since its discovery in neutral kaon decays in 1964 pQ. So far CP violation has 
been observed only in K and B hadron decays j2], results that are consistent with 
the Kobayashi-Maskawa formulation of a complex quark mixing matrix [3 , within the 
Standard Model (SM) of electroweak interactions. However, CP violation is not only 
important in hadron decays but it is also an inevitable ingredient for the generation 
of the observed baryon asymmetry of the Universe [3]. Yet, the observed CP violation 
within the SM is by far insufficient to explain it. This evidence is pointing to physics 
beyond the SM, whose contribution to CP violation in the K and B systems may be 
irrelevant but should be dominant for the Universe's baryon asymmetry. Thus, there 
are powerful motivations to study all possible CP odd observables and distinguish its 
underlying sources. Among those observables, a great deal of effort has been directed to 
the study of the electric dipole moment (EDM) of leptons, neutrons and neutral atoms. 
Both SM and various non-SM sources of CP violation have been considered (for recent 
reviews see e.g. Ref. jHl E])- These studies have been particularly stimulated by the 
expectation of great improvements (2 to 4 orders of magnitude) in the experimental 
sensitivities to EDMs in the next decade (for review see Ref. [B]). 

As it is well known, there are two sources of CP-violation within the SM: 
the complex phase of the Cabibbo-Kobayashi-Maskawa (CKM) matrix in the weak 
interaction sector, and the #-term in the strong interaction sector (Zj-^Hj. On the 
one hand, the complex phase of the CKM matrix is well established and provides a 
consistent explanation of the observed CP odd effects in hadron decays, but predicts 
an imperceptible contribution to the EDMs, far below the sensitivity of present or 
foreseeable experiments. Indeed, the CKM prediction for the neutron EDM, d n , ranges 
from 10~ 31 to 10~ 33 e-cm [TI] while its present experimental upper limit ^2] is 

|4| < 0.63 x KT 25 e-cm. (1) 

On the other hand, the presence of a #-term is still an open issue; it leads to a sizable 
neutron electric dipole moment [T?^ IT4] and may significantly contribute to atomic 
EDMs, while it is insignificant for CP violation in hadron decays. Indeed, the non- 
observation of EDM for the neutron imposes a very strict upper bound on the value 
of 6, of the order of 10~ 10 . This unnaturally small value of 6, which is otherwise not 
restricted by theory, is called the strong CP problem. One elegant solution was proposed 
by Peccei and Quinn (PQ) which makes the 9 parameter to vanish dynamically. 
Other important contributions to the atomic and neutron EDMs may arise from physics 
beyond the SM. 

In the calculations of atomic and neutron EDMs one faces the problem of translation 
of the CP violation mechanism formulated at the quark-gluon level, to the processes 
at the hadronic or atomic level. This translation must resort to hadronic and nuclear 
models and perform a careful treatment of the effects from the cloud of atomic electrons. 
The validity of the estimates then depends on the validity of these models. 
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The purpose of the present work is to use the PCQM model to calculate the 
neutron EDM and electric dipole form factor (EDFF) induced by the standard CP- 
odd QCD #-term as well as by generic CP-violating effective interactions arising from 
physics beyond the SM. The latter correspond to effective quark-photon and quark- 
gluon interactions induced by intrinsic electric and chromoelectric dipole moments of 
the constituent quarks and the Weinberg gluon operator at a more fundamental level. 
In the present paper we do not discuss concrete mechanisms for the generation of these 
CP-odd interactions, but concentrate on the hadronic structure aspects of the derivation 
of the neutron EDM in terms of these generic CP-violating parameters at low energies. 

The problem of the neutron EDM has been studied within different theoretical 
approaches: current algebra and chiral perturbation theory ^JQ^IEIEIj chiral quark 
models PSj-fS], lattice QCD [22l EHj, QCD sum rules [2U EEj, an approach based 
on solutions of Schwinger-Dyson and Bethe-Salpeter equations j2E], etc. We use the 
PCQM which describes baryons as bound states of three relativistic valence quarks 
confined in a static potential and supplemented by a cloud of pseudoscalar Goldstone 
bosons, as required by chiral symmetry. This model has already been successfully 
applied to the charge and magnetic form factors of baryons, ground state masses of 
baryons, the electromagnetic N-A transition, meson-nucleon sigma-terms and other 
baryon properties [2"T| . 

This paper is organized as follows. Sect. 121 contains a brief review of the PCQM. 
In Sect. 0] we apply the PCQM calculation of the EDFF and the EDM of the neutron 
induced by the individual EDMs of the quarks and their CP-violating couplings to the 
chiral fields. In Sections H] we extract constraints for the corresponding quark EDMs 
and CEDMs from the existing experimental limits on the neutron EDM and discuss the 
phenomenological implications of our results. 

2. The Perturbative Chiral Quark Model 

The basis of the perturbative chiral quark model (PCQM) |2Z| is an effective chiral 
Lagrangian describing the valence quarks of baryons as relativistic fermions moving in 
an external field (static potential) 

V cS (r) = S(r) + 7°V(r), with r = \x\, (2) 

which in the SU(3)-flavor version are supplemented by a cloud of Goldstone bosons 
(n,K,ri). Treating Goldstone fields as small fluctuations around the three-quark core, 
the linearized effective Lagrangian is written as: 

1 8 

C eS {x) = q(x)[i - S(r) - -y°V(r)]q{x) + - ^[^(x)] 2 

i=i 

+ £f (1 \x) + C x s B (x), (3) 
where we explicitly defined the linearized coupling of Goldstone bosons to quarks as: 

Cf r{1 \x) = -q(x) il ^S(r) q (x). (4) 
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The additional term C x sb in Eq. (J3J) contains the mass contributions both for quarks 
and mesons, which explicitly break chiral symmetry: 

£ x sb(x) = -q{x)Mq{x) - ^Tr[&(x)M] . (5) 

8 

Here, <E> = $jAj is the octet matrix of pseudoscalar mesons, F = 88 MeV is the 

i=l 

pion decay constant in the chiral limit, M. = diagjm, m,m s } is the mass matrix 
of current quarks (we restrict to the isospin symmetry limit m u = md = in) and 
B = — < 0|mm|0 > / F 2 is the quark condensate constant. We rely on the standard 
picture of chiral symmetry breaking and for the masses of pseudoscalar mesons we use 
the leading term in their chiral expansion (i.e. linear in the current quark mass): 

2 

Ml = 2rhB, M 2 K = (m + m s )B, M 2 = -(m + 2m s )B . 

In our analysis we use the following set of parameters: 

rh = 7 MeV, m s = 25m, B = M 2 +/2m = 1.4 GeV . (6) 

The meson masses satisfy the Gell-Mann-Oakes-Renner and Gell-Mann-Okubo relations. 
In addition, the linearized effective Lagrangian fulfills PCAC The properties of baryons, 
which are modeled as bound states of valence quarks surrounded by a meson cloud, are 
then derived using perturbation theory. At zeroth order, the unperturbed Lagrangian 
simply describes a nucleon as three relativistic valence quarks which are confined by 
an effective one-body static potential V e s(r) in the Dirac equation. We denote the 
unperturbed three-quark ground-state as \<fio), with the normalization (0 o |0o) = 1- We 
expand the quark field q in the basis of eigenstates generated by this potential as 

l( x ) = baUa (%) exp(-i£ a t) (7) 

a 

where the quark wave functions {u a } in orbits a are the solutions of the Dirac equation 
including the potential V e g(r). The expansion coefficients b a are the corresponding single 
quark annihilation operators. All calculations are performed at an order of accuracy 
0(l/F 2 ,rh,m s ). In the calculation of matrix elements, we project the quark diagrams 
on the respective baryon states. The baryon states are conventionally set up by the 
product of SU(6) spin-flavor and SU(3) C color wave functions, where the nonrelativistic 
single quark spin wave function is replaced by the relativistic solution u a (x) of the Dirac 
equation. 

In our description of baryons we use the effective potential of Eq. (j2J), which is given 
by a sum of a scalar potential S(r) providing confinement and the time component of a 
vector potential 7°V(r). Obviously, other possible Lorenz structures (e.g., pseudoscalar 
or axial) are excluded by symmetry principles. It is known from lattice simulations that 
a scalar potential should be a linearly rising one and the vector potential is thought 
to be responsible for short-range fluctuations of the gluon field configurations [2*%] . In 
our study we approximate Kfj(r) by a relativistic harmonic oscillator potential with a 
quadratic radial dependence 

S(r) = Mi + cir 2 , V(r) = M 2 + c 2 r 2 . (8) 
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The model potential defines unperturbed wave functions for the quarks, which are sub- 
sequently used to calculate baryon properties. This potential has no direct connection to 
the underlying physical picture and is thought to serve as an approximation of a realistic 
potential. Notice that this type of potential was extensively used in chiral potential 
models [29J- 31J. A positive feature of this potential is that most of the calculations can 
be done analytically. As was shown in Refs. and later on also checked in the 

PCQM [27 , this effective potential gives a reasonable description of baryon properties 
and can be treated as a phenomenological approximation of the more fundamental forces 
dictated by QCD. 

The use of a variational Gaussian ansatz for the effective potential (jHJ gives the 
following solution for the ground state (for the excited quark states we proceed by 
analogy) : 

u (x) = N exp ( \ ± J XsXfXc (9) 

where iV = [7r 3 / 2 i? 3 (l + 3p 2 /2)] _1 / 2 is a normalization constant, and Xs, Xfi Xc are the 
spin, flavor and color quark wave functions, respectively. The parameter p, setting the 
strength of the "small component" , can be related to the axial charge qa of the nucleon. 
In the leading order (3-quark-core) approximation, this relation is [2Z| 



5 / 2p 



2 



^-A l -TTW)- (10) 

The parameters of the effective potential V e g can also be expressed in terms of p and R: 

Mi = , Mo = On , Ci = Co = o • (11) 

2pR ' ° 2pR ' 2R 3 K J 

Here, £ is the single-quark ground-state energy. In our calculations we use the value 
(7a=1-25. Therefore, we have only one free parameter in the model, R. In our numerical 
study, R is varied in the region from 0.55 fm to 0.65 fm, which is set and constrained 
by nucleon phenomenology j22j. Such a variation of the parameter R slightly changes 
the physical quantities up to 5% 27 . In this paper we also test the sensitivity of the 
neutron EDM to a variation of R. 

The expectation value of an operator A is defined as 

(A)= B (4 >0 \f2- ] j d 4 Xl ... j a?x n T[C I {x 1 )...C I {x n )A]\<j )Q )* (12) 

n=l n ' J J 

where Cj is the full interaction Lagrangian which may contain both CP-even and CP- 
odd terms, as discussed below. The superscript "£>" in Eq. (|T2j) indicates that the 
matrix elements are projected on the respective baryon states and the subscript "c" 
refers to contributions from connected graphs only. 

For the evaluation of Eq. (fT2J) we apply Wick's theorem with the appropriate 
propagators for quarks and mesons. For the quark field we use a vacuum Feynman 
propagator for a fermion in a binding potential. In the calculation of meson-quark loops 
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we include only the ground state in the quark propagator, which leads to the following 
truncated form: 

iG q (x,y) = (0\T{q(x)q(y)}\0) -> 6(x - y )u (x)u (y}e-^ Xo ~^ . (13) 

Notice that in our previous papers [27 j we estimated explicitly the contribution of 
the low-lying excited quark states in the quark propagator to the physical quantities. 
Their contribution is about 10 to 15% with respect to the ground state contribution. 
Therefore, a restriction of the quark propagator to the ground states is a reasonable 
approximation. For completeness we also calculate the corrections to the neutron EDM 
due to the inclusion of excited quark states: the first p-states (lpi/2 and in the 

non-relativistic notation) and the second excited states (le/3/2, le/5/2 and 2si/ 2 ), i.e. we 
restrict to the low-lying excited states with energies smaller than the typical scale of 
A = 1 GeV of low-energy approaches. 

For the meson fields we use their free Feynman propagators: 

iA^x -y)= <0|T{$ l( *)$,(,)}|0> = J { C k )H %^:^ - ( 14 ) 

3. The neutron electric dipole Form Factor and Moment in the PCQM 

Here we apply the above described PCQM model to the calculation of the electric dipole 
form factor (EDFF) and moment (EDM) of the neutron. In this model, having the 
constituent quarks and pseudoscalar mesons as basic degrees of freedom, the effective 
chiral CP-odd Lagrangian at low energies can be written in the following generic form 

4pv(z) = ~\ fa) d q 0^ 75 q(x) F^(x) + 1 q(x) 75 h q e***^/^ q(x) 

= ~ 2 d q a ^ 7s q(x) F^ v (x) + % q(x) 75 h q q(x) 

- fa) {h q , H*)} fa) + 0($ 2 ) , (15) 

where F^ u is the electromagnetic field strength tensor, the symbol { , } denotes 
anticommutator, d q = diag(d u , dd, d s ) and h q = diag(/i u , hd, h s ) are diagonal matrices 
of the effective EDMs of the constituent quarks and their CP- violating couplings to the 
chiral fields, respectively. In Eq. (|15|) we expanded the quark-meson term in powers of 
&/F and kept the leading non-trivial term linear in $. It turns out that the first term 
in this expansion does not contribute to the EDFF at one loop because the diagrams 
involving this vertex do not contain the required spin-flip structure ■ q, the product 
of the neutron spin operator and of the 3-momentum of the photon q. Therefore, in 
the following we keep only the linear term in $ in Eq. fll5j) . 

We assume that the quark-meson Lagrangian (fTB^) originates from a more funda- 
mental level of quark and gluon CP-odd interactions, so that the effective parameters 
d q and h q can be related to those more fundamental CP-odd parameters. We discuss 
these questions in section |U 
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Now, we calculate the neutron EDFF, D n (Q 2 ), starting from the Lagrangian (|15|) . 
The EDFF is defined in the standard way through the neutron matrix element of the 
electromagnetic current as 

(n(p')\J»(0)\n(p)} = u n {p')W F£\Q 2 ) + -i- ^ q v F^{Q 2 ) (16) 

- a^ u 75 q v D n (Q 2 ) + (<y" q 2 - 2m N q») l5 A n (Q 2 )ju n (p), 

where, in addition to D n (Q 2 ), f£ ] {Q 2 ) and F^\q 2 ) are the well-known CP-even 
neutron electromagnetic form factors and A n (Q 2 ) is the neutron anapole moment form 
factor. The neutron EDM, d n , is defined as the value of the neutron EDFF at zero 
recoil, namely d n = D n (0). 

To guarantee electromagnetic gauge invariance in non-covariant approaches (see 
discussion in Refs. p?2"l l2Tj ) we have to work in the Breit frame, which defines the 
momenta of the transferred photon and the initial and final neutron states respectively 
as 

q=(0,q), p = (E,-q/2), p' = (E,q/2), (17) 

where E = m 2 N + q 2 /4 is the nucleon energy, the nucleon mass, and q 2 = —Q 2 = 
—q 2 the momentum transfer squared. We can now use Eq. ()16|) in the Breit frame and 
our model w.f. of Eq. (JHJ) to extract the form factor D n (Q 2 ) as the term in (J (0)) 
proportional to <Jat • q, namely 

E 
m N 

Here xn s and xjv , are the nucleon spin w.f. in the initial and final state, and 



U°(°)> = —xi,^N-qX Ns D n (Q 2 ) + .... 



(J°(0)> = (n(p')lf^ [ ^i-- - / d i x n T[Cr(x 1 )...C I (x n )J°(0))\n(p)) c , (18) 

n=l U ' J J 

where Lj is the full interaction Lagrangian and the subscript "c" refers to connected 
graphs only. In PCQM, the electromagnetic current operator 

J»(x) = Jg(x) + J» PV {x) (19) 

contains a CP-conserving electromagnetic part, given in terms of the charged pseudo- 
scalar meson fields: 

J£(ar) = e \ tt~ (x) id»ir + (x) + K~{x) id^K + (x) 1 + h.c. (20) 

and a CP-violating part involving the quarks: 



J c » PV (x) = i d v \q(x) d q 75 q(x)\ , (21) 

which is derived from the Lagrangian p5j1 . The interaction Lagrangian Li in Eq. (|18j) 
up to one loop order is the linearized strong interaction Lagrangian between quark and 
chiral fields, which consists of the CP-conserving interaction given in Eq. (j3J) and the 
CP- violating Lagrangian of Eq. (|15|): 

Cj{x) = -q( x )ijS^S(r)q(x) - ^ q(x) {h q , $(x)} q(x). (22) 
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Up to one loop, the neutron EDFF receives contributions from three types of effective 
diagrams, shown in Fig. 1: the "/i-terms" (Fig. la), which are proportional to the CP- 
violating couplings h q of Eq. (j22J), the "tree- level terms" (Fig. lb) and the "one-loop 
corrections" (Fig. lc), both due to the CP- violating current of Eq. (|2*Tj) and proportional 
to the couplings d q . The diagrams in Figs, la and lc encode the chiral corrections to the 
neutron EDM. The diagram in Fig. la has contributions from charged pion and kaon 
clouds, while the diagram in Fig. lc contains contributions from ir, K and //-meson 
clouds. 

The total contribution of Fig. la to the neutron EDFF is then composed by charged 
pion and Kaon clouds as: 

D%\Q 2 ) = D n' M (Q 2 ) > D [ n M (Q 2 ) = d^F^\Q 2 ) (23) 

where the normalization F^ h '^(Q) = 1 defines d\ h '^ as the contributions of the charged 
it and K clouds of Fig. la to the neutron EDM: 



j[h,<5>\ _ nih-Maw — \h-M 9tN n 9&n n f d 3 k F 7rNN (k 2 ) F nNN (k 2 ) 
d n -D n W -ec ^ J (27r)3 ^ 

and F^ h '^(Q 2 ), the form factor normalized to unity at zero recoil is given by: 



(24) 



j[h;$]/Q2\ _m N f d 3 k F wNN [(k + q) 2 ] F wNN (k' 



E J (2tt) 3 w$(k + q) w$(k) w^{k + q) + w$,{k) 

Here $ = n or K, w<s>(q ) = a/ M| + q 2 is the meson energy, and c^'*' is the SU(6) spin- 
flavor factor, which is c^ ;7r ^ = 1/2 for the pion-loop and c' h;i ^ = 1/5 for the kaon-loop 
diagram. 

The expression for d!n'^ in Eq. (}2*4*j) is written in terms of the known ttNN coupling 
which satisfies the Goldberger-Treiman relation 

9nNN = gA-^r, (26) 

(where g& = 1-25 is the axial nucleon charge) and the CP-violating couplings g<s>NN, 
which are: 

h u + ha _ h u + h s . . 

9*nn = 2F 7 , 9knn = 2F 7- I 2 ') 

Here 7 is the isovector-scalar two-quark condensate in the nucleon, given by: 

(N\qT 3 q\N) = ju N r 3 u N . (28) 

This factor 7 coincides with the so-called relativistic reduction factor [27J. In the PCQM, 
7 = 5/8 j2Zj- Finally, F nNN and F wNN are the normalized CP-conserving and CP- 
violating form factors that appear in PCQM |2Z|, arid which regularize the divergent 
loop integral: 

F nNN (k 2 ) = exp(-P J R 2 /4) [l + ™ (l - ^-)] , (29) 
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and 

K NN (k 2 ) = exp(-^ 2 /4) [l + k2 flgj!-\ ] ^ 

Now, it is worth checking if our result for the pion-cloud contribution to dn'^ in Eq. (|24j) 
is consistent with the model-independent prediction derived in Ref. for the leading 
term in the chiral expansion: 

d - = 6 4„*m N lQg A4- (31) 
To this end, we drop the normalized form factors F wNN and F nNN in Eq. (j^) by 
substituting F nNN = F nNN = 1, and analyze this equation using alternatively cutoff or 
dimensional regularizations. Both methods of regularization give the same result, which 
also coincide with Eq. (p!Tj) . Therefore, our approach is consistent with QCD in the local 
limit, when F n ^N = F n ^N = 1- The nontrivial form factors F n NN and F n NN provide an 
ultraviolet convergence for the EDM. Notice that Eq. (J3*T]) has been derived in ChPT 
when the strong CP- violating 7r iViV coupling constant g^NN is defined by the ^-term. 
This result is also valid for general CP- violating couplings. The difference is encoded in 
the redefinition of 

9nNN(0) = #^7 -> g^NN{h u ,h d ) = u ^ d j. (32) 

The leading contribution of the kaon-cloud diagram is also proportional to the chiral 
logarithm but contains a model-dependent coefficient c^ h ' K ^ = 1/5: 

flWCl = c [h;K] e 9«nn 9«nn lQg ^v (33) 
An^niN Mk 

We remark that the coefficient S- h ' ,K \ was calculated previously in Heavy Baryon Chiral 
Perturbation Theory (HBChPT) [17 , where it was expressed through the parameters 
of the chiral Lagrangian as: 

Jh-,K] _ F) — F bp — bp 

D + Fb F + b D [6) 
Here D and F are the axial- vector couplings, and bp, and bp are low-energy constants. 
Using the actual values \Tj\ D = 0.80 and F = 0.46 fixed from a fit of semileptonic 
hyperon decays, and b D = 0.079 GeV -1 , bp = —0.316 GeV -1 determined from the 
calculation of baryon masses and the 7T-N sigma-term up to fourth order in the chiral 
expansion, we deduce the prediction of HBChPT for d h ' K i = 0.45. This is more than a 
factor two larger than the prediction of our model. 

Now, let us examine the diagram Fig. lb. In this calculation we use the well-known 
quark model formula for the projection of the spin-flavor part of the Lagrangian (JT5j) 
between neutron states: 

<A 4 1 

(n t | 2_j ( d s); (° • 9 )i \ n t) = ~ z du ' ( 35 ) 

i=l 

where the sum is over the constituents, d q is a diagonal 2x2 flavor matrix of the quark 
EDMs, q = q/\q\ is the unit vector along the momentum q, and a are the Pauli matrices 
in spin space. 
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Then the contribution of the diagram in Fig. lb (a tree-level contribution) to the 
neutron EDFF is given by 

D^ ree] (Q 2 ) = - ^) Q + ^9a) FW(Q>) , (36) 



FM(Q 2 ) = ^exp(-Q 2 i? 2 /4) 



Q 2 R 2 1 - 3g A /5 
4 1 + 3W5 



rfid-Mtn^ - Jd;*] 9 *nn F [d]( 2\ f d3k p 
n W ' ~ 200 m 2 N W ' J (2vr)3 



$ = Tr,K,r). 



where F^(0) = 1. 

Finally, we present the loop correction diagram of Fig. lc. This one has three 
contributions: from 7r, K and r) loops, which we denote as: 

DW°°v\Q 2 ) = D^{Q 2 ) + D^ K \Q 2 ) + D^\Q 2 ) , (37) 

and where each meson loop contribution is given by: 

Here c [d ^ ] are SU(6) spin-flavor factors: c [d > n] = 7d u + 2d d , c [d '' K] = 6d s and c [d < v] = 
4d u /3 - d d /3. 

The total neutron EDFF, up to one loop in the chiral expansion, is then given by 
the tree level contributions of Fig. l.b and the loop contributions of Figs, l.a and l.c: 

D n (Q 2 ) = Df™\Q 2 ) + D^(Q 2 ) , 

Df ee \Q 2 ) = D^ tTee \Q 2 ) , (38) 

Octree] ^ = ^tree] ^ _ 1^ F M (g9) ^ 

Dj**] ( g2 )=aW (^4 + ^) ^ ] (Q 2 ), 

D^](g 2 ) = a^ 4 ^i(g 2 ), 
^ ] (Q 2 ) = « H (^-^) F [dl (Q 2 ), 



where 



(39) 



x 

where in the contributions of the diagram in Fig. la, for convenience, we use the 
dimensional parameter, A x = 4irF n ~ 1.2 GeV, the scale of spontaneously broken chiral 
symmetry, where F n = 92.4 MeV. The corresponding numerical coefficients a\ tree ^ = 7/8, 
a W = 3 x 10" 2 , = 6.3 x 10~ 3 , = 8.5 x 10~ 4 , &M = 2.92 and 6^ = 0.10 are 
obtained from the evaluation of Eqs. (|24|). ()36|) and (|37jl . 



The neutron electric dipole form factor in the perturbative chiral quark model 



11 



It is convenient to separate the contributions from the individual quarks as: 

D n (Q 2 )= D n,q{Q 2 ) = £ [A?d 9 ^(g 2 ) + Af^ g FfI(Q 2 )],(40) 

q=u,d,s q=u,d,s 



where we define h q = (e/A 2 ) h q and 
FuKQ 



, b^F^{Q 2 ) + b^F^ K \Q 2 ) 



6W + bW 

Ff\Q 2 ) = Fte*\(?) , (41) 
Ff](Q 2 )=F^(Q 2 ). 



The coefficients A q ^ and Ag' 1 ' are: 

A H = _ i a [*ree] _ ? M _ 4 W = 

3 9 3 

A W = 1 [«™.]_2 a W + I W = i.i 6 

d 3 9 3 

A [ f = - a [A1 = -0.006 , (42) 

AM = i&M + -6^ = 1.51, 
2 2 

A M = 1 6 W = L46 ) 

AM = ±&[*1 = 0.05. 
2 

The numerical results for the #-term contribution to the neutron EDFF (including the 
pion and kaon loops) are shown in Fig. 2. In Fig. 3 we plot the Q 2 -dependence of the d 
quark contributions to the neutron EDFF arising from its individual EDM, dj, as well 
as its CP-violating coupling to mesons, ha, contained in the first and second terms of 
Eq. pPjl. respectively. Since the values of dd and h q and their relative signs are not 
known, we display these two types of contributions separately, in terms of normalized 
form factors: F^ d (Q 2 ) = A^ F^(Q 2 ) for the individual d quark EDM contribution, 
and F^\(Q 2 ) = Ajj F^(Q 2 ) for the d quark contribution via its CP- violating coupling 
to mesons. 

The neutron EDM, d n , is then obtained by definition as the neutron EDFF at zero 
recoil. Using F^(0) = 1 and Ff ] (0) = 1, for q = u } d } s, the resulting EDM can also be 
expressed in terms of the individual quarks EDMs d q and CP- violating couplings h q , as: 
dn = £ [AWd i + A q h r hq ], (43) 

q=u,d,s 

with the coefficients A q and A q h ^ given in Eq. (|42p. 

Another CP-odd parameter which may have important implications for CP- 
violation in atoms is the electron-neutron Schiff moment 5", related to the neutron 
EDM form factor D n (Q 2 ) by [33] 

'dD n (Q 2 ) 



S' 



dQ 2 



(44) 

Q 2 =0 
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From the expression (|38|) - (j42|) we obtain: 

gl _ g> [h] _|_ gl [d] 



(45) 



where 
S' [h 



e ' 

& = TT,K 



2m N 



d 3 k F nNN (k 2 ) F nNN (k' 



(2tt)' 



w%(k 



x 



5 k' 



+ 



Ml F^ NN {k 2 ) 



m 



N 



IP 

w%{k) F wNN (k 2 ) 3 /•; 



■NN 



(k' 



nNN 



(k' 



is the contribution of the diagram in Fig. la and 



S'ld] 



X 



3 



10 



9a 



$=ir, K,rj 



fjA 



5m 



N 



2 1 + 3g A /5 



q=u,d,s 



Aid] A 



[d;$] 9-kNN 


1 



200 



d 3 £; ^ 2 F 2 NN (k 2 ) 



N 



W W%{k 

+ - R 2 9A 
2 1 + 3^/5 _ 



is the contribution of the diagrams in Figs, lb and lc. Here F^ NN and F" NN are the 
first- and second-order derivatives of the F wNN form factor with respect to k 2 . 

Again, as for the case of the EDM, we can show that our Schiff moment is consistent 
with the ChPT results in leading order of the chiral expansion. As in the case of 
the neutron EDM, we drop the normalized form factors F^nn and F^nn in Eq. (J2HJ), 
substituting F nNN = F tNN = 1 and, therefore, F^ NN = F" NN = 0. Then we analyse 
the leading term in the chiral expansion (which is ultraviolet convergent): 



S' 



e ' 



9ttNN9<$>NN 

2m N 



A straightforward calculation of the integral 
S'=e £ 



d 3 k 1 
gives 



5 k 2 

6 w%(k 



e ^ c 



[h;<£>] 9tvNN g<S>NN 

48 7T 2 m N M 2 



(46) 



(47) 



In the case when we restrict Eq. (jjTj) to the contribution of the pion-cloud, our result 
coincides with the leading order result of ChPT |33| I34j : 

9-kNN 9-kNN 



S' 



(4f 



48 vr 2 m N M 2 ' 

Let us now present our results for the electron-neutron Schiff moment S' in terms of the 
partial contributions: 

(a) Diagram in Fig. la: partial 7r and K meson loop contributions 



gl [dy 



16.22 GeV~ 



x 



-h u + -h d 



S f[h;K] = q >21 Ge y- 



\hu 



2 * 



(b) Diagram in Fig. lb: tree-level result 

S ,[d;tree] = g^g Ge y-2 



X 



— dd — —do, 
3 3 



(49) 



(50) 
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(c) Diagram in Fig. lc: partial tt, K and rj meson loop contributions 
s >ldi*] = _ .31 GeV~ 2 x 



7 2 " 
9 u + 9 d 



S '[d;K] = _ Q 06 Ge y-2 x d ^ ^ 
S '[d;v) = _ .01 GeV~ 2 x 



— do, — — dr! 

3 3 



Next we summarize our result for the Schiff moment in terms of the individual quarks 
EDMs d q and CP-violating couplings h q , as: 

S'= [P ld] ^ ] d q + (3 q h] A^h q }, (52) 

q=u,d,s 

with the coefficients A q ^ and given in Eq. (}4*2*|) . Here, for convenience we introduced 
the /9-factors: (3^ = — ^[^'(0) and (5^ = —Fq'(0). Their numerical values are: 

= 10.06 GeV" 2 , # = 5.44 GeV" 2 , (if = 5.55 GeV~ 2 , pf ] = 2.10 GeV~ 2 . (53) 
4. Limits on the quark EDMs and their phenomenological implications 



Here, using Eq. ()43|) we extract the limits on the quark EDMs, d q , and CP-violating 
couplings h q , from the existing experimental upper limit on the neutron EDM, shown in 
Eq. (Q). Assuming the absence of accidental cancellations between the different terms 
in Eq. ()43)1 we obtained the corresponding limits shown in the Table 1. 

In order to discuss possible phenomenological implications of these limits and 
compare them with the results of other approaches, it would be desirable to express 
the individual quark EDMs in terms of more fundamental CP violating parameters 
specifying the origin of CP violation. 

The effective CP-odd Lagrangian in terms of the quark and gluon fields up to 
operators of dimension six has the following standard form [3*3]-|39j: 

C ^ {x) = J^ tl{dfluGn ~ 2~ g ~ d * a " u 75 F » u q ~ l q d ? °" v 75 G ^ T ° q 

-\c w f abc G% G b « G% , (54) 

where G a is the gluon stress tensor, G^ u = ^^vopG 011 is its dual tensor, and T a and 
jabc are foe SU(3) generators and structure constants, respectively. In this equation 
the first term represents the SM QCD 0-term, while the last three terms are the non- 
renormalizable effective operators induced by physics beyond the SM. The second and 
third terms are the dimension-five electric and chromoelectric dipole quark operators, 
respectively, and the last term is the dimension-six Weinberg operator [33]. The quark 
EDMs, = diag(dy , , df ), and quark chromoelectric dipole moments (CEDM), 
d^ = diag(d^, d%, d^ ), form diagonal matrices in flavor space. The operators of the 
above Lagrangian can be induced by physics beyond the SM after integrating out the 
heavy degrees of freedom. 
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The contribution of the QCD #-term can be related to the parameters of the 
Lagrangian (JT5|) on the basis of the following standard reasoning. A chiral U(l) 
transformation in flavor space allows one to remove the gluonic #-term from the 
Lagrangian ([54)1 and pass it as a CP-violating complex phase to the quark mass 
operators. For a small value of 6, the CP- violating term becomes (for details see 
Refs. [CH EJ HOI): 

C s cpv=idrh Q{ x )lsQ{x), (55) 

q=u,d,s 

where q(x) denotes a triplet in flavor space, and so this term is a flavor-5'?7(3) singlet. 
The mass coefficient is: 

m u m d m s 

m = , (5o) 

m u m d + m u m s + m d m s 

which would vanish if any of the quarks were massless. From this term one can construct 
an effective chiral Lagrangian by introducing the chiral field e* 75 *^ and expanding it in 
powers of $/F: 

£ S CPV = iQmq~l5 exp f i^% ) q = iOmq^q - Orhq^q + 0(<l 2 ) . (57) 



For the reasons we gave after Eq. (|15|). the leading contribution to the EDFF comes 
from the term linear in $ in the expansion of Eq. (|57|) . Comparing the Lagrangians in 
Eqs. (fT3j) and (JSTj) we obtain the contribution of the (9-term to the parameter h q which 
is, naturally, flavor-blind: 

h q = hf M + 9ml 3x3 , (58) 

where I3X3 is the unit matrix in flavor space. Here we have singled out the part 
hq SM which does not contain the 9 contribution but is due to the other terms in the 
Lagrangian (|54j) . 

In order to relate the parameters d q and h? of the Lagrangian (|T3jl with the 
parameters d^, d^ and Cw of the Lagrangian we apply a simplified approach 
relying on the so-called naive dimensional analysis [23 El H2] • I n this way one obtains 
the following order-of-magnitude relations: 

d, =< + + ^ C »'- ( 59 > 

h r M =^, + ^c w . (60) 

We assume all the operators in Eq. (J54|) to be normalized at an energy scale of about 
1 GeV, where the perturbative quark-gluon picture is still reasonable and, on the other 
hand, occurs an overlap with the chiral symmetry breaking scale A x = AttF^ ~ 1.2 
GeV, which corresponds to the normalization scale of the Lagrangian (|15|) . Thus, all 
the parameters in Eqs. ()59|) and (jUU|) are at the same scale A x ~ 1.2 GeV. In such a case, 
it is also natural to identify the parameter A with A x . In a specific model beyond-SM, 
the parameters of the effective Lagrangian (p>4*j) could be calculable in terms of the model 
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parameters at some high energy scale. These parameters then must be QCD-evolved 
down to the hadronic scale A x corresponding to the scale of the neutron EDM. The 
renormalization-group relations of the parameters d^, d^ , Cw in Eqs. (|5T?|) and (|HU|) 
with their values calculated at a high energy scale can be found in Refs. jHSl EH] • 

Now, using Eqs. (JHHj). ([671]) in Eq. (JI3J), we can express the neutron EDM in 

terms of the more "fundamental" parameters 9, dg, d^ and Cw of the Lagrangian (jSljl 
as 

= ( z f d f + ) + Zw A x ^ + (^«x 1(T 16 e • cm) . (61) 

The numerical values of the ^-coefficients are given in Table 2, where we also present for 
comparison the values of these coefficients derived in some other approaches: QCD sum 
rules (23 El|, parton quark model |4~K) . SU(6) model jlHj, MIT bag model [IB], current 
algebra [Tlj, effective chiral approach JJj], HBChPT chiral bag model [TP] , cloudy 
bag model [20], chiral quark- meson model [21] and extensions of SM jUj-jSH]- 

Notice that our results for the coefficient z e is dominated by the pion-loop diagram. 
The kaon-cloud contribution is smaller that the pion contribution by a factor ~ 1/28: 

z e = z* e +zf = 1.42, (62) 
4 = 1.37, zf = 0.05. 

Using Eqs. ([58)1 . ([BT?]) . (jfiTj)) in Eq. (JS2J), we can write for the electron- neutron Schiff 
moment the following expression: 

S' = (sf tZf + ) + s w A x C w + (6 s e x 10~ 16 e ■ cm) . (63) 

where numerical values of the s-coefficients are 

sf = - 4.86 GeV" 2 , sf = 17.85 GeV~ 2 , sf = -0.09 GeV~ 2 , 

s° = 0.41 GeV" 2 , s c d = 1.62 GeV~ 2 , sf = 0.004 GeV" 2 , (64) 

s w = 2.04 GeV -2 , s =7.72, 

and where the separate contributions of pion and kaon meson-loop diagram in Fig. la 
to the coupling se are s$ = 7.62 and sf = 0.1. 

Now, using Eq. ([51)1. we can extract upper limits for the parameters (9, d 9 , h q , df, 
d^ and CV from the current experimental bound on neutron EDM in Eq. ([!)). These 
limits are shown in the Table 1. In their derivation we assumed the dominance of 
each term at a time, thus assuming the absence of significant cancellations between the 
different terms contributing to the neutron EDM. 

We note that, in comparison with the conventional valence quark model which 
predicts for the neutron EDM 
4 1 

d n = ^d d - -d u , (65) 

the PCQM predicts a small but non-negligible contribution from the strange quark 
EDM, which appears due to the K-meson loop diagram in Figs, la and lc. In Table 2 
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we show the predictions for the neutron EDM of some other approaches. With the 
values of the corresponding ^-coefficients it is straightforward to derive upper limits on 
the CP-violating parameters of the Lagrangian (|54jl within the cited approaches. A 
general summary is that our limits are less stringent. Notice that all the approaches 
cited in Table 2, including the PCQM, were successful in describing various properties of 
baryons, and yet they disagree in the CP- violating sector, as they predict different values 
of the neutron EDM. This a manifestation of a strong hadronic model dependence in the 
estimates of the neutron EDM. In view of that, the existing EDM limits on fundamental 
CP violating parameters could only be treated as a hint for model building rather than 
as stringent limits. 

5. Summary and Conclusions 

In this work we applied the perturbative chiral quark model to the calculation of the 
electric dipole form factor and electric dipole moment of the neutron. We parameterized 
a generic effect of CP violation in terms of the effective quark electric dipole moments 
and strong CP-violating quark couplings to the chiral fields, which could have different 
sources at some more fundamental level. We considered the relations between these 
effective parameters of the quark-meson Lagrangian with more fundamental parameters 
of the CP-odd quark-gluon Lagrangian. The latter includes the SM #-term and the 
operators which originate from physics beyond the SM, such as the chromoelectric 
and electric quark dipole moment operators, as well as the Weinberg term. From the 
existing limits on the neutron EDM we extracted the limits on the individual effective 
quark EDMs and on the aforementioned more fundamental CP violating parameters, 
adopting the approach based on the so called naive dimensional analysis. A general 
conclusion we made from the comparison of these results with other approaches is that 
our limits are significantly less stringent. This could be treated as a manifestation of a 
strong dependence of the predictions for the neutron EDM on the details of the hadronic 
structure models. 
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Table 1. Upper limits on the CP- violating parameters of the Lagrangians in 
Eqs. (|15|l . (|54|l and CPVCs h q ; the high-energy constants d^ , d^ and from the current 
experimental bound on the neutron EDM in Eq. JQ). All numbers are given in units 
of 10~ 25 e-cm. 



Parameter 


u 


d 


s 


\dg\ 


2.0 


0.5 


105.0 


\hq\ 


0.4 


0.4 


12.6 


\df\ 


1.3 


0.4 


68.6 




6.4 


2.9 


174.7 


\9\ 


0.4 


\C W A x \ 


2.0 



Table 2. Theoretical estimates of z-coefficients in Eq. (|61ll in various approaches. 



Parameter 


Approaches 


This work 


z E 


- (0.35 ± 0.15) gU ; -0.78 g3] ; -1/3 gB] 


-0.48 




(1.4 ± 0.6) gU ; 1.14 gS] ; 4/3 gH] 


1.77 




-0.35 gS] 


-0.01 


z c 

u 


0.48 [IH! ; (0.165 ± 0.075) [HJ ; 0.05 |Hj ; 


0.10 




4/9 e ~ 0.135 gE] 




% 


0.39 [H] ; (0.33 ± 0.15) gl] ; 0.05 gTj ; 


0.21 




8/9 e ~ 0.270 gH] 




*° 


0.08 [H] ; 0.05 glj ; 0.03 48J ; 0.008 49 


0.004 


z w 


~ 0.085 [SH! ; ~ 0.017 glj ; ~ 0.01 50J 


0.094 


z 9 


2.7 [13]; 3.6 [H|; 3.3 [Ej; 6.7 [IZj; 3.0 [HI!; 


1.42 




1.4 20j 1.17 EH; 2.4 [23] 





i* i i * i r i 

(a) (b) (c) 



Figure 1. The contributions to the neutron EDM: "/i-term" diagram (a), "tree-level 
ci-term" diagram (b) and "one-loop d-term" diagram (c). The black filled squares 
correspond to the CP-violating vertices of the Lagrangian ljT5)) or the current if^Tl) . 
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Figure 2. The contributions from pion and kaon loops to the neutron EDM form 
factor induced by strong CP-violating #-term, where Dn\Q 2 ) = Dn l ^(Q 2 ) at the limit 
h q = 9 fh I 3x3 . 




0.1 0.2 0.3 

Q 2 [GeV 2 ] 



0.4 



Figure 3. The <i-quark contribution to the neutron EDM form factor via its intrinsic 
EDM (solid line) and CP-odd coupling to mesons (dashed line). The plotted form 
factors are: F^(Q 2 ) = A l * ] F^(Q 2 ) (solid line) and F^Q 2 ) = A^ 1 Ff 1 (Q 2 ) 
(dashed line). 



